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ABSTRACT 

We present 37 GHz data obtained at Metsahovi Radio Observatory in 2001 December - 2005 
April for a large sample of BL Lacertae objects. We also report the mean variability indices and 
radio spectral indices in frequency intervals 5-37 GHz and 37 - 90 GHz. Approximately 34 % 
of the sample was detected at 37 GHz, 136 BL Lacertae objects in all. A large majority of the 
detected sources were low-energy BL Lacs (LBLs). The variability index values of the sample 
were diverse, the mean fractional variability of the sample being AS'2=0.31. The spectral indices 
also varied widely, but the average radio spectrum of the sample sources is flat. Our observations 
show that many of the high-energy BL Lacs (HBL), which are usually considered radio-quiet, 
can at times be detected at 37 GHz. 

Subject headings: BL Lacertae objects: general - galaxies: active - radio continuum: galaxies 



1. INTRODUCTION 

BL Lacertae objects (BL Lacs, BLOs) demon- 
strate the most violent behaviour known among 
active galactic nuclei (AGN). Their basic proper- 
ties include high optical polarisation, rapid and 
strong variability, and a nea r-featur eless spectrum 
( Stein et al. 1 119761: iKollgaar d 1994; iJannuzi et al.l 
1994, for an extensive review on AGN, see 



Urrv &: Padovanil (jl995[ )). The selection criteria 
in many surveys also include a flat radio spectrum 
(a < 0.5, S oc v-") (jPadovani fc Giommil [l995a . 
and the references therein). 

The spectral energy distributions (SEDs) of BL 
Lacertae objects are thought to consist of two com- 
ponents; a synchrotron component at lower fre- 
quencies and an inverse Compton (IC) component 
at higher frequencies. The locations of the compo- 
nents in the log - log vFi^ plane varies much, and 
this has been used as a classification basis for BL 
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Lacs ( Padovani fc Giommi' 'l995b; I Giommi et al 



19951 : ISambruna et al.Hl996 : Nieppola et al.ll20o'a 



hereafter Paper I). For low-energy peaked BL Lacs 
(LBLs) the peak frequency of the synchrotron 
component is at radio and IR-frequencies, for in- 
termediate BL Lacs (IBLs) in the optical and UV- 
frequencies, and for high-energy peaked BL Lacs 
(HBLs) in the X-ray band. To further explore the 
physics behind this spectral sequence, it is impor- 
tant to know which are the properties that set 
LBLs and HBLs apart. 

Radio monitoring provides insight into the vari- 
ability characteristics and long-term behaviour of 
AGN. Metsahovi is one of the few observatories 
which target the high radio frequencies, and the 
only one observing at 37 GHz. BL Lacs have 
been observed in Metsahovi as part of the ongo- 
ing Metsahovi - Tuorla Observatory collaboration 
project. There are currently 398 BL Lacs in the 
Metsahovi AGN sample. We wanted to study this 
large sample of BL Lacs at 37 GHz because many 
of them had not previously been studied at high 
radio frequencies. Our aim was to measure the 
flux densities of the various subsamples of BLOs 
in order to get a full understanding of the high 
radio frequency behaviour of these objects, all the 
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way from the radio-selected BL Lacs (RBLs) to 
the X-ray selected BL Lacs (XBLs). HBLs and 
XBLs have not been actively studied in the high 
radio frequency domain. Even though XBLs seem 
to be significantly weaker at radio wavelengths 
than RBLs, the re is l ittle e vidence of radio-silent 
BLOs. istocke erZI (|l99d ) suggest that they do 
not exist, and even though there have been some 
recent studies on pos sibly completely rad i o-silen t 
BL Lac candidates (jLondish et all I2OO2I |2004) . 
the number of such BLOs is probably low. Be- 
cause the only high radio frequency s tudy o n XBL s 
is a very small sample studied by 
and only the low frequency radio flux densities 
of these sources are known in general, it is pos- 
sible that some of the rarely studied BLOs ex- 
hibit "surprising" spectral shapes in the high radio 
frequency domain, such as the extreme in v erted- 
spectrum sources (Tornikoski et al. 20001 



2001 



Torniainen et al. 20051 ) 



Our studies also support the work of Extra- 
galactic foreground sources Working Group of the 
Planck satellite0, particularly that of the Low Fre- 
quency Instrument (LFI) Consortium. In prepa- 
ration for the Planck mission, one of our tasks is 
to estimate the number of extragalactic sources 
detectable at the Planck frequencies. The exis- 
tence of unforeseen bright sources could seriously 
affect the primary task of Planck, i.e., the map- 
ping of the cosmic microwave background (CMB). 
The theoretical aspects of blazar contamination 
in CMB maps have been thoroughly studi ed by 
Giommi fc Colafrancesco ( 20041) : Giommi et al 



(20061): Toffolatti et al 



(199S 



Because of the 
lack of actual measurements at high radio fre- 
quencies, many of the source statistics have been 
based on extrapolated low-frequency data. There- 
fore we find it important to investigate, also in 
practice, if there are source populations or sub- 
populations that are brighter in the high radio 
frequencies than earlier assumed, especially such 
that exhibit significant variability and can at times 
make a significant contribution to the contamina- 
tion of the CMB maps. It is vital to understand 
the high radio frequency behaviour of the various 
BLOs to see how many of them, at least some of 
the time, can also be detected by Planck. 

The main objective of this paper is to publish 
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the data we have collected for almost three and a 
half years. We have also calculated some descrip- 
tive parameters for the sample. In §2 we introduce 
the sample and describe the observations, in §3 we 
state the detection rates by class. In §4 and §5 
we report the values of the variability indices and 
spectral indices we found for the sample, and in §6 
say a few words about the variability timescales. 
We finish with a discussion and conclusions in §7 
and §8. All statistical tests have been performed 
with the Unistat 5.0 software. 

2. THE SAMPLE AND OBSERVATIONS 

The Metsahovi BL Lac sample consists of 398 
objects, 382 of which are from the Veron-Cetty 
& Veron BL Lac catalogue (Veron-Cetty & Veron 
2000, hereafter VCV2000), and 17 objects from 
the literature. IVCV200fll consists of sources from 
several different surveys with a wide range of se- 
lection criteria. The Metsahovi sample is one of 
the largest BL Lac samples ever studied and of- 
fers a relatively unrestricted view of this group 
of AGNs. The declination range covered by the 
Metsahovi BL Lacs is from -11° to 86°. 

Most of the sample source s were classified 
as LBL/IBL/HBL in iPaper 11 based on their 
synchrotron peak frequencies, Vpeak, calculated 
from a parabolic fit to their SEDs. The class 
boundaries were log Vpeak 

< 14.5 for LBLs, 14.5 < 
logVpeak < 16.5 for IBLs and logVp^ak > 16.5 for 
HBLs. There are 98 LBLs, 96 IBLs and 110 HBLs 
in the sample. 94 sources remained unclassified 
on account of poor spectral fits. The full source 
sample is listed in Table [TJ For most sources, we 
use the ob ject naming convention adopted from 
VCV20001 The subclass of each source is also 
listed. 

Metsahovi radio telescope is a radome enclosed 
antenna with a diameter of 13.7 metres. It is sit- 
uated in Kirkkonummi, Finland, at 60 m above 
the sea level. It has been used for AGN variabil- 
ity studies at 22, 37 and 87 GHz. The observa- 
tions for this project were carried out using the 37 
GHz receiver, which is a dual horn, Dicke-switched 
receiver with a HEMT preamplifier operated at 
room temperature. The observations are ON-ON 
observations, alternating the source and the sky in 
each feed horn. A typical integration time for one 
fiux density data point is 1200-1600 seconds, and 
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Table 1 

The Metsahovi BL Lac sample. 



Source name Right Ascension Declination class'' 





(h) 


(min) 


(s) 


(deg) 


(arcmin) 


(arcsec) 




NRAO 5 





6 


13.9 


-6 


23 


36 


LBL 


RX J0007.9+4711 





7 


59.9 


+47 


12 


7 


IBL 


MS 0011.7+0837 





14 


19.7 


+8 


54 


4 


HBL 


RXS J0018.4+2947 





18 


27.8 


+29 


47 


32 




PKS 0017+200 





19 


37.9 


+20 


21 


46 


LBL 



Note. — Table [T]is published in its entirety in the electronic edition of the Astronomical 
Journal. A portion is shown here for guidance regarding its form and content. 

''as determined in lNieppola et al.l l|2006l ). 



Table 2 

37 GHz FLUX DENSITIES FOR THE METSAHOVI BL LAC SAMPLE. 



Source name Date of the observation (J2000) Flux density Error 







(year) 


(month) 


(day) 


(h) 


(min) 


(Julian date) 


(Jy) 


(Jy) 


NRAO 5 




2002 


5 


24 


6 


40 


2452418.778 


1.71 


0.17 


NRAO 5 




2003 


3 


30 


9 


20 


2452728.889 


2.94 


0.25 


NRAO 5 




2005 


2 


9 


13 


10 


2453411.049 


2.69 


0.16 


RX J0007.9- 


h4711 


2002 


3 


17 


6 


28 


2452350.769 


<0.36 




RX J0007.9- 


h4711 


2003 


3 


16 


5 


56 


2452714.747 


<0.2 





Note. — Table [2] is published in its entirety in the electronic edition of the Astronomical Journal. A portion 
is shown here for guidance regarding its form and content. The error is listed for S/N > A detections only; for 
non-detections we have calculated S/N > 4 upper limits. 
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the detection limit under optimal weather condi- 
tions is approximately 0.2 Jy. As a primary flux 
calibrator we use DR21 and as secondary calibra- 
tors 3C84 and 3C274. The error values include the 
contribution from the measurement rms and the 
uncertainty of the absolute calibration. For more 
details about the Metsahovi observing syst em and 



data reduction see Terasranta et al. ( 1998h . 



The observations reported in this paper were 
performed between years 2001.95 and 2005.27. All 
flux densities and upper limits are listed in Ta- 
ble [2l The number of observations, iVobs, for each 
source varies greatly. The most observed sources 
include BL Lac, OJ 287, Mrk 421, AO 0235+164 
and S5 0716-f 714. They aU have more than 100 
data points, S5 0716-1-714 topping the list with 
355 measurements. Such a high number of mea- 
surements is exceptional because an overwhelm- 
ing majority of the sample, 96%, have A'obs < 10. 
The source S5 0716-1-714 was the target of an in- 
tens ive multifrequency ca mpaign in 2003 Novem- 



ber ( Ostorero et al. 20061 ). Before the start of the 



actual campaign it exhibited extremely active be- 
haviour in the radio domain and was monitored 
very frequently especially in late 2003 and 2004. 

The observation interval for each source varies 
also; for sources with iVobs > 5 the average interval 
is 0.34 years, or approximately 124 days. For ob- 
jects observed 5 times or less, the average interval 
is of the order of 13 months. 

Some of the best-known, radio-bright sources in 
the sample have been observed at Metsahovi also 
prior to 2001 as par t of the extragalactic sources 



moni t oring project ( Terasranta et al. 1992 . 19981 



2004 120051 ). Many of them have frequently- 
sampled flux curves dating back to the beginning 
of 1980's. If these data are included, 21 sources 
in all have more than 100 measurements. We will 
discuss the flux curves and variability timescales 
of this limited BL Lac sample using the full data 
set in a forthcoming paper (Nieppola et al. 2007, 
in preparation). 

3. DETECTION RATES 

Table [3] lists the detection rates for the whole 
sample and the BL Lac classes separately. The de- 
tection limit is S/N > 4, typically ranging from 0.2 
to 0.6 Jy depending on weather conditions. Ap- 
proximately one third of the sample has been de- 
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Fig. 1. — The distribution of the average fluxes of 
the sample sources. 

tected at S/N > 4 at 37 GHz, although all of them 
have been measured at least once. Most of the de- 
tected sources are LBLs, and only one tenth are 
HBLs. This is due to the larger average radio lumi- 
nosities of LBLs compared to HBLs. Luminosity 
correl ations of BL Lac classes are discussed in de- 
tail in lPaper j . The detected HBLs include objects 
which were at first deemed unlikely to be detected. 
Typically, there are a few non-detections and one 
detection. This suggests that there is variability 
in these sources as well, but normally we cannot 
detect it due to the faintness of the source. Fig. [1] 
depicts the distribution of the average fluxes of the 
sample. It emphasizes the difficulties of detecting 
the sources; roughly half of the detected sources 
have average fluxes under 0.5 Jy. 

4. VARIABILITY AMPLITUDES 

To examine the variability amplitude of the 
sample we calculated 3 variability indices: 

— {Sniax S rfiiiyi) I S jnin (1) 
— {Sjnax Sfjiin) I {,Srnax + Sjyiiji) (2) 
— Smax ~ Sfnin- (3) 

The flrst two describe the fractional variability, 
whereas the third is simply the absolute variability 
of the source. AS*! gives explicitly the proportion 
of the flux change to the minimum flux. Thus 
AS'i = l means that the source has increased its 
flux level to twice its minimum flux density during 
our observing period. 
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Table 3 

Detection rates for the Metsahovi sample. 



Class 


number of 
detected sources 


% of class 
detected 


% of all 
detected sources 


number of 
undetected sources 


% of class 
undetected 


% of all 
undetected sources 


all sources 


136 


34.2 


100 


262 


65.8 


100 


LBL 


75 


76.5 


55.1 


23 


23.5 


8.8 


IBL 


35 


36.5 


25.7 


61 


63.5 


23.3 


HBL 


16 


14.6 


11.8 


94 


85.5 


35.9 



Note. — The detection limit is S/N > 4. All detected sources have at least one S/N > 4 measurement. For all sources, A^obs > 1. 



Only S/N > 4 detections were used in calculat- 
ing the indices. The bare minimum of two detec- 
tions needed to estimate the 37 GHz variability 
were available for 87 sources in total, of which 6 
were HBLs, 17 IBLs, 58 LBLs, and 6 were un- 
classified. Fig. O shows the distributions of the 
indices. The average and median variabilities for 
the whole sample and BLO classes separately are 
listed in Table [H 

In the distributions of ASi and AS2 (Fig. ^ 
upper and middle panels) most of the objects can 
be found in the lower end of the range, followed 
by a high-variability tail. A clear exception from 
the general distribution is S5 0716+714, whose 
AS*! is anomalously high, 17.47, and AS'2==0.9. 
For this object an unprecedented radio outburst 
was recorded during our observing period (see also 
§6). The peaks of the distributions are AS*! =0.5- 
1.0 and AS2 =0.2-0.3. Looking at the lower panel 
of Fig. [21 it is clear that although the fractional 
variability of the sample sources is relatively high, 
their absolute variability is not. Even out of the 
detected sources, a large majority is very faint and 
the flux levels increase less than 0.5 Jy. Notewor- 
thy are the three sources in the far end of the scale, 
exhibiting flux changes of more than 5 Jy. These 
objects are S5 0716-^714, OJ 287 and BL Lac. 

As we see from the average results in Table |4j 
LBLs seem to exhibit slightly stronger variability 
than HBLs. The variability of IBLs is approxi- 
mately the same as that of HBLs. However, there 
is a significant correlation between the number of 
detections, A'det and variability (Fig. [3]), according 
to Spearman Rank Correlation test. Radio-bright 
LBLs are observed more often and therefore they 
get a larger A^det- The more measurements there 
are of a source, the more probable it is that it is ob- 



served both during a quiet and a flaring state, thus 
increasing its detected variability. To eliminate 
this selection effect, we calculated revised variabil- 
ity parameters for a sample in which A^det < 9 for 
all sources. For those 30 sources that had more 
than 9 detections, 8 representative points were se- 
lected, with as regular time intervals as possible. 
We chose Adet=8 as the limit of number of data 
points because below that point the correlation 
between A^det and AS2 in Fig. [3] disappears. The 
results for the Afdet-hmited sample are also listed 
in Tabled! 

Limiting the value of A'det affected the two ex- 
tremes, LBLs and HBLs, the most, while the vari- 
ability parameters of IBLs remained almost the 
same (Table 01). This is because there are several 
LBLs with extremely high A'det and therefore very 
high variability, the most extreme object being S5 
0716-1-714 for which Afdet=353. Two examples of 
the six HBLs in the sample are the radio- luminous 
Mrk 421 and Mrk 501 (Afdet=74 and A^dct=38, re- 
spectively). Because of the small number of HBLs, 
limiting the A'det-values of these two objects had a 
big impact on the average variability index of the 
sample. 

We looked for a possible correlation between 
AS2 (Fig. m and the s ynchrotr on peak frequency 
i^peak (as calculated in jPaper J ) using the Spear- 
man Rank Correlation Test. There was no sig- 
nificant correlation (P =0.25). The lack of cor- 
relation between A5'2 and the synchrotron peak 
frequency holds also for the A'det -limited sample. 
This further emphasizes the continuity and uni- 
formity of the BLO population found in Paper I. 
Usually RBLs are said to be more variable than 
XBLs, but our results so far give no sign of such 
a distinction. 
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Table 4 

Variability indices for the whole Metsahovi sample and BL Lac classes separately. 



Class 


ASi 


A^2 


AS3 




A^2 


A^3 


Ndct 


Ndct 




average 


average 


averag 


e median 


median 


median 


average 


median 


all data 


all sources 


L19 


0.31 


0.69 


0.84 


0.30 


0.39 


18.05 


5.00 


LBL 


L36 


0.33 


0.87 


0.93 


0.32 


0.41 


21.22 


6.00 


IBL 


0.84 


0.26 


0.32 


0.50 


0.20 


0.24 


10.00 


4.00 


HBL 


0.78 


0.26 


0.32 


0.83 


0.30 


0.34 


20.67 


3.00 


iVdot < 9 


all sources 


0.85 


0.27 


0.55 


0.72 


0.26 


0.38 


5.41 


5.00 


LBL 


0.93 


0.29 


0.68 


0.83 


0.29 


0.39 


5.74 


6.00 


IBL 


0.70 


0.23 


0.28 


0.50 


0.20 


0.20 


4.88 


4.00 


HBL 


0.52 


0.20 


0.24 


0.50 


0.20 


0.27 


4.67 


3.00 



Note. — Indices are as described in SJH and A'dot is the number of S/N > 4 detections for 
one source. In the upper half all data is included, while in the lower half data is limited to 
^dct < 9. See text for details. 



5. SPECTRAL INDICES 

We calculated spectral indices for frequency in- 
tervals 5-37 GHz and 37-90 GHz to investigate 
the spectral behaviour of BL Lacs below and above 
37 GHz. The spectral index is defined as oc v". 
Only S/N > 4 detections were used. The results 
by class are shown in Table [5] The 37 - 90 GHz 
results for HBLs are omitted beacuse there were 
only 2 HBLs with both 37 and 90 GHz data. 
The distributions of the indices can be found in 
Fig. [SI The 5 GHz data are taken from the As- 
trophysi cal Gatalogues S upport System (CAT^, 
IVcrkho danov et aljll997h maintained by the Spe- 
cial Astrophysical Observatory, Russia. The 90 
GHz data are taken both from CATS and from 
our observations with the Swedish-ESO Submil- 
limetre Telescope (SEST) between 1987 and 2003. 

In Table [H] and Fig. [5] there are 3 different in- 
dices, aquiet, ct flare and aave- They are calculated 
from the minimum flux densities, maximum flux 
densities, and average flux densities, respectively. 
In both frequency intervals the range of indices 
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was quite large. The maxima of average indices 
were 2.23 and 1.09, and minima -0.56 and -1.47 
in 5-37 GHz and 37-90 GHz, respectively. As we 
see from Fig. [SJ the majority of the sources have 
flat (a >-0.5) spectra in both intervals. The 5 - 
37 GHz distribution has a tail towards the higher 
indices. There are 18 objects with aave > 1- Typ- 
ically, these sources are X-ray selected BLOs with 
little radio data. In 37 GHz in particular, 13 
sources have only one S/N >4 datapoint. It is 
plausible that the source has been in a flaring state 
at the time of the observation, which would result 
in an inverted spectrum. The 37 - 90 GHz distri- 
bution is less skewed but broader. S5 0716+714 
has a conspicuously steep spectrum; its aave —- 
1.47. The steepness is exaggerated by the non- 
simultaneity of the data; the 37 GHz average flux 
is strongly influenced by the high flux levels of the 
2003 flare, whereas the 90 GHz data is from a qui- 
escent state at an earlier epoch. For the whole 
sample, most of the data taken at 37 GHz were 
from different observing epochs than the 5 GHz or 
90 GHz data. Thus we possess no information of 
the true shape of the continuum spectra during the 
various activity states. For a small subset of our 
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sample we have managed to obtain truly simul- 
taneous multifrequency spectra between 1 and 37 
GHz. These results will be discussed in Tornikoski 
et al. (in preparation). The non-simultaneous 
spectral indices calculated in this paper only give 
us a crude estimate of the spectral behaviour. 

There are no major differences in the spec- 
tral indices between the BLO classes among the 
detected sources. For HBLs the 5-37 GHz in- 
dices are slightly more inverted than for LBLs and 
IBLs. Spearman Rank Correlation test found no 
correlation between the spectral indices and the 
synchrotron peak frequency Vpeak, except for a 
marginal positive correlation for afiare and aave 
in 5-37 GHz. 

Generally, there is a correlation between spec- 
tral index and variability, steep spectrum sources 
being on ave rage less var i able than flat spec - 
trum sources ( Edelsonlll987t Valtaoia et al. 1992[ ). 
However, among the flat spectrum sources, and es- 
pecially among BL Lacs, no significant correlation 
has been detected. Our study makes no exception: 
the 37 GHz variability depends on neither of the 
spectral indices according to the Spearman Rank 
Correlation test. 

6. TIMESCALES AND FLUX CURVES 

While all detected BL Lacs seem to exhibit 
some degree of variability, the timescales are di- 
verse. We calculated the slopes of the flux density 
changes for each source. S5 0716-1-714 has clearly 
the fastest flux density variations. It has exhibited 
documented intraday variability (IDV) behaviour, 
and has also been a subject of a very intensive 
monitoring campaign in 2003. Its flux density in- 
creased by 1.06 Jy in 40 minutes and 1.12 Jy in 
less than 2 hours, and decreased by 0.87 Jy in less 
than 2 hours during the 2003 flare. The multifre- 
quency flare and the IDV during the c ampa ign are 



discussed in detail in lOstorero et al. 

AO 0235-1-164 has the second fastest variations: 
its flux density increased by approximately 0.4 Jy 
in less than an hour in 2004 August, after drop- 
ping from 2.38 to 1.76 Jy in 40 minutes 16 hours 
earlier. Also OJ 287 is worth mentioning with its 
flux density increase of 0.8 Jy in approximately 2 
hours. 

There were 5 sources, which had a well-sampled 
flare during our observing period: AO 0235+164, 



S5 0716-^714, OJ 287, 1308-^326, and BL Lac, all 
of them LBLs. No two flares were alike; their type 
ranged from powerful and short (S5 0716-1-714, 
2003 July - 2003 December, peak flux density 6.28 
Jy) to broad and weak (1308-1-326, 2002 January 

- 2004 October, peak flux density 3.28 Jy). OJ 
287 had two different outbursts. The later one 
in 2003 August - 2005 April reached the highest 
flux density value recorded in our observations, 
7.8 Jy. The flare of BL Lac was still ongoing in 
2005 April. The flux density had risen steadily 
since 2004 February, and had reached 6.6 Jy in 
2005 March. The mean duration of a flare is 17.6 
months. 

Most of the sources in our sample have very few 
data points and thus no flux curve to speak of. As 
mentioned in S}3l for many of these faint sources 

- mainly HBLs and XBLs - our data consist of 
detections (S/N> 4) and non-detections (S/N< 
4). This is partly due to the varying observing 
conditions, with the detection limit ranging from 
0.2 to over 0.5 Jy, but it can also, at least partially, 
be caused by variability. There are examples of 
sources, expected to be faint and usually below 
the detection limit, suddenly measuring over 0.5 
Jy, such as MS 0737.9-1-7441 and GB 1011-1-496. 
An especially good example of such behaviour is 
RXS JlllO.6-1-7133. It has one detection at 0.74 
Jy and 6 non-detections. 

7. DISCUSSION 

In this paper we publish the data of the first 3.5 
years of our ongoing BL Lacertae object observing 
project. We have also calculated some variability 
parameters and spectral indices for the sample. 

The sample consists of 398 BL Lacs. During 
the 3.5 years of monitoring we made at least one 
significant {S/N > 4) detection of shghtly more 
than one third of the sources at 37 GHz, with 
the detection limit being of the order of 0.2-0.5 Jy. 
The largest subpopulation that was detected was 
LBLs (77%). 37% of IBLs and 15% of HBLs were 
detected. Even though the percentage of the de- 
tected BLOs is smaller in the high-energy peaked 
population than in the low-energy peaked one, it 
is interesting to see that also among them there 
is a large number of objects that were detected 
at 37 GHz. This has a significant implication 
for the extragalactic foreground point source mod- 
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Table 5 

Spectral indices by class for frequency intervals 5 - 37 GHz and 37 - 90 GHz. 



Class 


O^quiet 


aflart 




O^quiet 


flare 


O^ave 




average 


averag 


e average 


median 


median 


median 


5-37 GHz 


all sources 


0.33 


0.17 


0.25 


0.16 


0.00 


0.05 


LBL 


0.30 


0.15 


0.23 


0.15 


-0.02 


0.01 


IBL 


0.28 


0.07 


0.16 


0.13 


-0.02 


0.02 


HBL 


0.67 


0.51 


0.60 


0.68 


0.55 


0.51 


37-90 GHz 


all sources 


-0.10 


0.15 


0.00 


-0.12 


0.22 


-0.02 


LBL 


-0.13 


0.21 


0.01 


-0.10 


0.36 


-0.01 


IBL 


-0.03 


0.23 


0.09 


-0.18 


0.22 


0.13 



Note. — 37-90 GHz indices of HBLs are omitted due to small sample 
size. The indices have been calculated with non-simultaneous archival 
data. 



elling of the Planck satellite. Because the Planck 
detection limit will be comparable to that of the 
Metsahovi telescope, approximately one third the 
BLO population could also be detected by Planck. 
This should include at least some of the high- 
energy peaked BLOs that are often excluded from 
both high-frequency radio observing campaigns as 
well as many of the models for CMB contribu - 
tion by blazars. Giommi &: Colafrancescol 1 20041 ) 
suggested that blazar flux variability at millime- 
tre wavelengths may be a significant contributor 
to the CMB map contamination. Our findings of 
variability within the BLO sample, including clear 
signs of variability for sources in the HBL popula- 
tion, confirm this. 

In order to have a reasonable estimate of how 
probable it is to observe the faint BLOs above 
the detection limit, we calculated the ratio of de- 
tections versus non-detections. For the HBLs for 
which there is at least one detection in our data, 
but excluding the two densely monitored and usu- 
ally detected objects Mrk 421 and Mrk 501, the 
percentage of detections among all good quality 
data points was 19%. If we add the X-ray selected 
BLOs from the IBL and LBL subclasses, the per- 
centage is 22%. Thus the probability to detect an 



HBL or XBL at a random observing epoch is of 
the order of 20%. This either means that their 
flux density remains on the verge of our weather- 
dependent detection limit, or that they are vari- 
able, flaring sources. 

Our earlier studies on Southern fi at-spectrum 
radio sources ( Tornikoski et al. 2000f ) show that 



many radio-bright AGNs spend much less time in 
an active state than in a quiescent or intermediate 
state. Our lon g term monitoring proj ect of radio- 
bright AGNs (jTerasranta et al.ll2005l . Hovatta et 
al. 2006, in preparation) shows that 3.5 years of 
data is a rather short time to make any conclu- 
sive statements about such variability behaviour 
as typical peak fluxes of flares or peak-to-peak 
timescales, because many of the dramatic changes 
in the flux density occur on long timescales. Thus 
it is important to consider the possible effects of 
long term variability when interpreting results of a 
few years' monitoring campaign. While we can say 
that the BLOs detected during this campaign can 
be bright in the mm-domain, we cannot rule out 
the possibility that the ones that now remained 
undetected could be much brighter if observed in a 
flaring state. This is an important issue when con- 
sidering the contamination of CMB maps caused 
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by extragalactic foreground sources, the data used 
for the modeUing are one to few-epoch historical 
data. We will discuss the long term variability of 
a subsample of these BLOs in more detail in Niep- 
pola et al. (2006, in preparation). 

The variability indices AS 2 of our sample range 
from 0.04 to 0.9. The brightest sources are ob- 
served more often whereas with the faint ones we 
make do with just a couple of measurements. This 
creates a selection effect and leads to biased vari- 
ability behav iour, as explained in f|4l For instance, 
in [ciaramella et al.l (|2004) the 37 GHz variability 
index for LBLs was AS'2=0.74 ±0.27. Our average 
value was AiS'2=0.34. Our value is lower because it 
includes many rarely observed faint LBLs. There- 
fore it is very problematic to assign any variability 
value to a source, much less a source class, as it is 
always case-specific. Considering the differences in 
the results of the complete sample and the iVdot - 
limited one, we find that the variability index AS2 
changes little compared to ASi and AS3. This is 
because, as the difference between Smax and Smin 
increases, AS2 approaches unity. Meanwhile, A5i 
grows without a limit and AS3 approaches Smax- 
Therefore AS2 seems to be least sensitive to the 
effect of a varying number of measurements, and 
as such probably the most convenient description 
of variability. 

When we compare the variability characteris- 
tics of the BLO classes - LBLs, IBLs and HBLs - 
no definite differences can be found. This is unex- 
pected. Radio-selected BL Lacs (RBLs) are gener- 
ally thought to be on average more variable than 
X-ray selected BL Lacs (XBLs). Our results indi- 
cate that this is not necessarily true for LBLs and 
HBLs. However, we have to keep in mind that 
we can only study the radio-bright HBLs, whose 
properties one could expect to be similar to those 
of LBLs. Over 80 % of HBLs still remain unde- 
tected at 37 GHz. 

The 5 - 37 GHz and 37 - 90 GHz spectral in- 
dices of the sample are diverse. In the lower fre- 
quency band the spectra are mostly flat, as ex- 
pected, but in the higher band the distribution 
broadens to include both rising and falling spec- 
tra. This indicates that although BLOs are flat 
spectrum sources at low frequencies, in the high 
radio frequencies their spectral behaviour is not 
necessarily uniform. 

The high-index tail of the 5-37 GHz spectral in- 



dex distribution includes mostly XBLs with only 
one detection at 37 GHz. This could be thought 
of as further evidence of the flaring nature of the 
faint BLO population. If we assume that in qui- 
escence, these XBLs have similar spectra to those 
of the vast majority of BLOs, their quiet flux level 
would be roughly of the order of 50 mJy or less. 
Thus, to be detected at Metsahovi, their flux den- 
sities must increase four- fold at the very least. We 
believe that active monitoring of the faint BLO 
population would result in more such detections. 
Also, using a more sensitive observing system to 
study a limited sample would be fruitful. 

8. CONCLUSIONS 

We have observed a sample of 398 BL Lacertae 
objects in Metsahovi Radio Observatory from 2001 
December to 2005 April. We obtained a S/N > 4 
detection for 34% of the sources. The detected 
fraction of LBLs was 77%, of IBLs 37% and of 
HBLs 15%. These figures are concordant with 
th e typical radio-faintn ess of HBLs, as observed 



Nieppola et al. ( 20061 ). However, we stress that 



HBLs cannot be considered radio-silent; this is 
verified by the serendipitous detections we ob- 
tained for some of them. 

We gratefully acknowledge the funding from the 
Academy of Finland for our Metsahovi and SEST 
observing projects (project numbers 205969, 
46341, and 51436). We acknowledge EC funding 
under contract HPRN-CT-2002-00321 (ENIGMA). 
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(jVerkhodanov et al.l 119971 ) of the Special Astro- 
physical Observatory. 
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Fig. 2. — The distributions of the variability 
indices AS'i, AS2 and AS3. The object S5 
0716+714 (A5i =17.47) has been omitted from 
the upper panel for the sake of clarity. 



CO 

< 



0.6 - • 

0.5 - . 

0.4 -. ':\ 

3 1; 

0.2 -'if.- 

0.1 'f ■ 



S5 0716+714 
AO 0235+164 



20 40 



60 80 



100 120 140 



Fig. 3. — Variability AS'2 plotted against the 
number of detections, A^det- Data points for 
S5 0716+714 (iVdct=353, A52==0.9) and AO 
0235+164 (iVdet=205, AS'2=0.44) have been ex- 
cluded from the figure for the sake of clarity. 
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Fig. 4. — Variability AS'2 plotted against the syn- 
chrotron peak frequency Vpeak ■ 
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Fig. 5. — Distributions of the quiet, flare and average spectral indices for the sample in the frequency 
intervals 5 - 37 GHz and 37 - 90 GHz. 
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